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SUMMARY 


Pressure-distribution expressions and stability derivatives have 
been derived by use of linear theory for zero-end-plate triangular 
vertical tails with subsonic leading edges performing rolling, yawing, 
and constant- lateral- acceleration motions. Corresponding results for 
the sideslip motion, most of which have been previously reported, are 
also included herein. 

Consideration is given to the effect of end plates on the forces 
acting on the vertical tail. Stability-derivative formulas for a 
vertical tail in the presence of a complete end plate obtained from 
wing results are also presented, together with a suggested approxima- 
tion for partial-end-plate effects. 

The aerodynamic damping of the lateral oscillation in yaw is 
approximated to the first order in frequency from the damping of the 
yawing and constant-lateral-acceleration motions. Illustrative varia- 
tions of the stability derivatives for the special case of the half- 
delta tail for all the motions considered are included. 


INTRODUCTION 


Information available at present that pertains to the aerodynamic 
forces acting on various tail arrangements is, in many instances, insuf- 
ficient to allow the accurate prediction of the lateral dynamic behavior 
of aircraft traveling at supersonic speeds. Theoretical results now 
available are, for the most part, concerned with tail configurations 
either in a rolling or in a sideslip attitude (see refs. 1 to 11). 
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For the sideslip motion, the effects of Mach number and aspect 
ratio on the aerodynamic loads of a number of tail configurations with 
both one and two planes of cross-sectional symmetry have already been 
investigated extensively. The same effects on tail arrangements in a 
rolling motion have also received considerable attention but it has, in 
the main, been directed toward tails with two planes of symmetry such 
as cruciform arrangements. Additional theoretical analysis devoted to 
the evaluation of the Mach number and aspect-ratio effects on the forces 
and moments acting on tail systems in roll with one plane of cross- 
sectional symmetry is required. 

Tail arrangements performing a steady yawing motion or a constant- 
lateral-acceleration motion have received little attention to date in 
the literature. Yet the forces and moments produced by these motions 
are by no means negligible, and some indication of their magnitudes is 
necessary, particularly at supersonic speeds, in order to evaluate their 
relative importance on lateral stability. 

The primary puiq)ose of this paper is to provide the pressure- 
distribution expressions and corresponding stability derivatives for 
isolated triangular vertical tails with subsonic leading edges per- 
forming yawing, rolling, and constant- lateral-acceleration motions. 

Some of these results in turn are used to approximate to the first 
order in frequency the damping of the vertical tail oscillating in 
yaw (one degree of freedom). 

A secondary objective, in view of the geometric nonplanar charac- 
teristics of tail arrangements, is to give consideration to the estima- 
tion of the mutual aerodynamic interference that exists between the 
vertical and horizontal tails. In this connection the stability deri- 
vatives for the vertical tail in the presence of a complete end plate 
have been included. 

For completeness, results for the no-end-plate and complete-end- 
plate vertical tails in a sideslip motion obtained from references 9 
and 11 are also presented. 


SYMBOLS 


The positive directions of the forces, moments, velocities, and 
angles are shown in figure 1. 

x,y,z coordinates of field point 


0 


z 

x 
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Xi^Zi coordinates of doublet 


a 


£1 


^o^^o 

t 

u 


V 

w 

V 
a 

P 

r 


a 

a 

3 

3 

p 


<1 

AP 


B = 



K 


distances origin is displaced relative to tail apex 
time 

incremental velocity in x-direction 

y-component of velocity 

z-component of velocity 

free-stream velocity 

speed of sound 

rolling angular velocity 

yawing angular velocity 

angle of attack 

rate of change of a with time 
sideslip aingle 

rate of change of 3 with time 
fluid density 

1 p 

free-stream dynamic pressure^ 

pressure difference between opposite sides of a surface 


1 

constant determining degree of homogeneity of quasi- 
conical velocity field 


0 

0D 


velocity-potential function 

potential of supersonic doublet distribution 
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potential of line of doublets 


A(x,z) 

f(cr) 

Sv 


steady-state potential corresponding to unit angle of 
sideslip 

steady-state potential corresponding to unit yawing 
velocity about z-axis 

doublet- strength f\inction 

line-doublet-distribution function 

span of vertical tail 

vertical- tail area 


e 

C 

R 


“"V _ 


2C 


1 - R 


aspect ratio of vertical tail, 

apex angle of tail 
tangent of apex angle 

ratio of slope of leading edge of tail to slope of 

2BC 


trailing edge of tail, 1 - 

BA-ii 


k = 


1 - \/l - b2c2 


BC 


E' (BC) 


K' (BC) 


complete elliptic integral of second kind vith modulus 

/l - b2c 2, r ^ ^1 - (l - B^c^jsin^n dn 

^0 


complete elliptic integral of first kind with modulus 

il - b 2 c 2 , r'^ ^ 

° /l - (l - B2c2)sin2n 


complete elliptic integral of second kind with modulus 

/ pJ( / 2 

l/l - k^, / |/l - (l - k^Jsin^n dn 

Jo 


E’(k) 
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K'(k) 


complete elliptic integral of first kind with modulus 


° /i - (1 - 


k^jsin^n 


m,n,m',n' 


constants 


N = NBC 
4 

G(BC) = 


7 = 


infinitesimally small quantity 
1 - 

(1 - 2b2c2)e'(BC) + B^C^k'CbCJ 

X - B^zff 


/l - B^cr^ /x^ - B^(y^ + z^j 
1 - B^ae 

I = lim 7 = 

)/l - B2a2/l - b202 


Cl 


rolling-moment coefficient^ 


■'ll 


yaving-moment coefficient, 


Rolling moment 
q^Syby 

Yaving moment 
cqSyby 


side-force coefficient, 


Side force 

q.Sy 


/Sc- 


3->0 


Cl. = 


'be- 




Cl. = 


V /p — 40 
\ V /r-40 
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Subscripts ; 

p rolling condition 

r yawing condition 

p sideslip condition 

p constant- lateral-acceleration condition 

0 zero-end-plate configuration 

C coraplete-end-plate configuration 

P partial-end-plate configuration 


SCOPE 


Derived in this paper are expressions for the surface pressure 
distributions on isolated trianguleir vertical tails performing yawing^ 
rolling, and constant- lateral-acceleration motions. These pressure 
expressions have in turn been used to ceilculate the stability deriva- 
tives associated with the side force, yawing moment, and rolling moment 
due to constant yawing ^^r)^ constant rolling 

Cnp^ ^<3. constant lateral acceleration 


Also presented are the aerodynamic coefficients Cy^, and 

Cjp obtained from the sideslip pressiire-distribution expression given 

in reference 11. Stability derivatives presented for the vertical 
tail of figure 1 mounted on a complete end plate are Cyp, Cnp, 


Cy^, 




Cy*, and 


Cn^. The equation of for a triangular- 

'p P 


vertical- tail — end-plate combination may be derived from the pressure 
distributions given in reference 9* This equation has not been deter- 
mined and only curves of Cjp taken from reference 9 fo^^ a vertical- 

tail — complete-end-plate combination with trailing edges perpendicular 
to the root chord (see fig. 2 ) are presented. 


The stability derivatives presented herein are valid within the 
limits of linear theory for a range of Mach numbers for which the leading 
edge is subsonic and the trailing edge supersonic. Variations of all 
the stability derivatives with the parameter BC are presented for a 
triangular vertical tail with zero trailing-edge sweep (half -delta tailj. 
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ANALYSIS 


The axes system used in the analysis and the positive directions 
of the forces, moments, velocities, and angles are shown in figure 1. 

The positive directions of the forces and moments in the analysis sys- 
tem have been fixed to conform to the positive directions of the forces 
and moments in the stability axes system (see fig. 5(a)). The stability- 
derivative expressions are derived in the body of the report with respect 
to the analysis system whose origin is at the apex of the tail. Transfer 
formulas are presented, however, which allow these derivative expressions 
to be determined for a system of axes whose origin is displaced longi- 
tudinally and vertically with respect to the apex of the vertical tail. 
(See fig. 5(b)-) 

The ensuing analysis is based on linearized theory and the results 
are restricted to vertical tails of zero camber with surfaces of van- 
ishingly small thickness. These conditions implicitly stipulate that 
the results are valid only for small angles of sideslip, small rates 
of change of the sideslip angle with time (p motion), and low rates of 
rolling and yawing. 

Throughout the analysis, when the potential of a surface in the 
xz-plane is referred to, it is considered to be the surface potential 
on the side of the vertical tail whose outward normal is in the posi- 
tive y-direction; that is , 


0 = 0(x,O'*',z) 

'and for a surface in the xy -plane, 

0 = 0(x,y,O'^) 

Pressure differences between opposite sides of a surface are formulated 
as follows: 

For a surface in the xy-plane, 

AP = P(x,y,0'*') - P(x,y,0“) 
and for a surface in the xz-plane, 

AP = P(x,0+,z) - P(x,0",z) 
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Determination of Pressure-Distribution Expressions 


for Yawing and Rolling Motions 


A method for solving supersonic- flow boundary- value problems 
governed by the classical, linearized, partial-differential equation 


has been developed in reference 12 and an application to rolling and 
pitching triangular wings is given in reference 15. This method allows 
the prediction of the disturbance-potential function 0, and hence the 
pressure distribution, for planar lifting surfaces. The analysis given 
in reference 15 is briefly summarized herein and is applied to the 
determination of the pressure distributions and associated forces and 
moments acting on a triangular vertical tall surface (fig. l) performing 
rolling and yawing motions. (Yawing in the xz-plane is analogous to 
pitching in the xy-plane . ) 

The determination of the form of the velocity potential .- As is 
well-known, the potentials of both the supersonic soiirce and the super- 
sonic doublet and their distributions represent solutions of equation (l 
For the determination of the potentials and pressure distributions of 
lifting surfaces of the type considered herein, that is, for lifting 
surfaces with subsonic leading edges , it is well-known that a distribu- 
tion of doublets that unlauely satisfies the prescribed bo\mdary condi- 
tions must be determined. These boundary conditions on the vertical 
tail for the motions to be considered herein are as follows: 

On the rolling vertical tail. 



( 1 ) 


Sy^ 


V = pz = xp ^ = xp0 


( 2 ) 


and on the yawing vertical tail. 


V = -rx 


( 5 ) 


In addition, 
of the tail: 


the following relations must be valid on the surfaces 
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For the rolling motion. 


< 5 ). 


ae 


p = p 


( 4 ) 


and for the yawing motion, 


^2/v 




P = 0 


ae‘ 


( 5 ) 


ae 


= 0 


( 6 ) 


aeP 


= 0 


( 7 ) 


The potential in space produced by a distribution of doublets, for 
example, in the xz-plane, with the doublet axes normal to the plane is 


0p(x,y,z) = 

ay 


-A(x-L,Zj^)dx2^ dzj_ 


/(x - X 3 _)^ - B^(z - Zp)^ - S2y2 


(8) 


where the area S is the region of the xz-plane intercepted hy the 
forecone from the field point (x^y^z). 

The potential on the surface carrying the doublet distribution is 
given by 


0p(x,Z)y^O 


lim 

y -^±0 


- ff- 

-A(x 2 ^,Z]_)dx]_ dz]_ 

ay 1 

/(x - Xp)^ - B^(z - zp)^ - s2y2 
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As stated in reference 15; this surface potential is directly propor- 
tional to the doublet-strength function a(x,z); that is, 

0D(x;z)y=±o ^ *rtA(x,z) (9) 


The surface-pressure velocity 

u(x,z)y=±o = 


u(x,z) 

d0jj(x,z) 


therefore is 
y=*0 _ 0 a(x,z) 

— Xjt 

Sx 


and the linearized lifting-pressure coefficient 

^ V 


may be written as 


^ ^ 4jt ^A(x,z) 

<1 V 3x 


( 10 ) 


( 11 ) 


( 12 ) 


The problem to be considered in this paper is one in which the 
sidewash on the surface is prescribed (see eqs. (2) and (5)) and the 
surface velocity potential has to be determined. The doublet- strength 
function A(x,z) then is an unknown and the determination of this 
quantity requires in general the solution of an integral equation. In 
some cases the general form of the surface-potential function A(x,z) 
is known or can be obtained by inverting an integral equation. The 
problem then resolves simply into an evaluation of the arbitrary con- 
stants of the general solution by making use of the prescribed boundary 
conditions . 

Brown and Adams in their analysis of triangular wings with subsonic 
leading edges (ref. 15) were able to determine the function A(x,z) for 
these wings undergoing various motions by utilizing the concept that the 
conical properties of the produced flow gave rise to potentials and 
pressures in the crossflow planes that were similar in form to the 
potentials and pressures acting on flat finite segments in a two- 
dimensional flow; these segments correspond to a section of the wing in 
any crossflow plane. This remarkable connection between linearized 
supersonic conical flow and incompressible two-dimensional flow is dis- 
cussed by Busemann in reference 1^4-. 
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A more general and rigorous approach to obtain the doublet- strength 
function may be formulated from an analysis presented in a later paper 
by Lomax and Heaslet (ref. I 5 ) dealing also with conical and the so-called 
quasi-conical problems. In this analysis a general surface-pressure- 
coefficient expression has bpen determined for planar lifting s\irfaces 
with prescribed boundary conditions of the form 


V ~ x«g(|) 


This expression is 



K+1 




i=0 /(C - e )(0 - Cq) 


(15) 


(14) 


where bi are constants, ® ~ determined by the boundary- 

condition equation (eq. ( 13 )), and C and Cq are the tangents of the 
apex angles of the two panels of the lifting surface. When Cp = 
the lifting surface is syinmetrical about the common root chord of the 
two panels^ and when Cp ^ the lifting surface is asymmetrical about 
this chord. From equation (12)^ which relates the function A(x^z) to 
the pressure coefficient^ and equation ( 9 ) the fom of a(x^z) or^ 
synonymously^ the form of the surface potential, may be obtained by a 
simple integration. It should be mentioned at this point that refer- 
ence 15 presents a method for deriving the arbitrary constants bj_ in 
the pressure coefficient (eq. (l^l-)). This method is related to that of 
reference 13 which concerns itself with obtaining the arbitrary constants 
in the velocity potential. 

By application of equation (l^) to the boundary problem of the 
vertical tall sketched in figure 1 (Cq = O) and by noting from the pre- 
scribed boundary conditions (see eqs. (2) and (3)) that k = 1, the 
pressure coefficient for both the yawing and rolling motions is 

^ = X ^0 + + bge^ 

^ ^ /(C - 0)0 

The constant bQ in this expression must be set equal to zero in order 
to satisfy the condition that along the streamwise edge the pressure 
must be zero. 
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The velocity potential on the vertical tail surface is easily 
obtainable from the pressure expression by the formula 




and has been found to be 


0=nx2fg) ( 16 ) 

where 

f(|) = f(0) = (»© + Nc)0e(c - 0) (it) 

The arbitrary constants in the so-called distribution function f(0) 
are, in terms of b^^ and b 2 ^ 


M = 


V ^ W\ 

^3c2 C J 


2itB 3 c2 

By relating equation (l 6 ) to equation (9)^ the doublet- strength func- 
tion A(x,z) is seen to be 

A(x,z) = x2f^|^ ( 18 ) 

A comparison of the potential of equations (l6) and ( 17 ) and the 
potential obtained for the slender, rolling, vertical tail reported in 
reference 5 shows, as expected, that both are of the same form. 

Evaluation of the constajits M and N .- The constants M and N 
in the expressions for the velocity potential given by equations ( 16 ) 
and ( 17 ) are still to be determined. As indicated previously, the 
expression for the pressure coefficient, and hence the velocity poten- 
tial, can be determined completely through an application of the pro- 
cedures developed in reference I 5 ; however, many of the integrations 
and Integrating procedures required in the method in reference 15 were 
a!Lready known to the author at the inception of this project and, for 
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this reason, the analysis herein to determine the constants M and N 
closely parallels the procedures discussed in reference 13 . 

The determination of the constants M and N depends upon satis- 
fying the boundary conditions associated with the vertical tail for the 
rolling and yawing motions . These bovindary conditions are given by 
equations (2) to (7). The needed expressions for the prescribed veloc- 
ities and their derivatives with respect to 0 in terms of the dis- 
tribution function f(cr) are derived in appendix A. 

For the rolling motion the constants M sind N may be obtained 
by replacing f(c) by its equivalent (eq. (17)) in the equations given 

in appendix A for v/x and and then applying the boundary 

conditions given by equations (2) and (4). When the integrations have 
been performed, the resulting equations may be solved simultaneously 
for Mp and Np. The yawing constants are obtained in a like manner 

with equations (3) and (6) replacing equations (2) and (4), respectively. 


In the calculation of the quantities v/x and 


any value 


^(v/x) 

ae 

of 9 may be considered. It is advantageous for integration purposes 
to let this value of 0 be zero. However, since one of the limits of 
integration is zero and since in several of the integrands a singular 
point exists at 9 = o = 0, the integrations in which these singularities 
occur must be performed for 0 arbitrary and then 0 is set equal to 
zero. 


Substituting equation (17) into equation (A4) gives, for 0 equal 
to zero. 



-3(MBa + NBC)\/Ba(BC - Ba)tanh"ll/l - ^ 2(MBa + NBC)/Ba(BC - Ba) 


(1 - b2o2) 


5/2 


d(B0) + 


(1 - b2ct2) 


lim 

0—^0 




lim ^ (MBa + NBC)|^a(BC - Bg)(l - B^oB)^/! - B^e^ 

B 


(1 - B 2 a 2 )^(Ba - B 6 )^ 


^ pBC (MBa + NBC)i/Ba(BC - Ba)(l - B^oB )^|/l - b202 


1 r 

® ^B(e+ii) 


(1 - b2c2)^(B0 - B0)^ 


d(Ba) - 


2(MB0 + NBC)\/Be(BC - B0 ) - 3^2 






( 19 ) 
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Carrying out the integrations in equation (l9) yields 

|K'(k) [ik^M + Nk2(l + 


V _ 
X 


b/i + k2(l - k2)' 


where 


and 


E'(k)[^(l + k2) - N(l - 4k^ + k^)]| 
1 - \Jl - b 2 c 2 


k = 


BC 


( 20 ) 


( 21 ) 


N = NBC = N — 

1 + k2 

These integrations were accomplished with the aid of the tables in 
references l6 and I7 and are discussed in appendix B. 

Substituting the distribution function into equation (A5) results 
in, for 9 approaching zero, 


^(v/x) _ 

pc 

3Bo(MBo + NBC)/bo(BC - Ba)tanh"^/l - B2q2 3 Bo(MBo + NBC)vfea(BC - Bo) 

de 

L 

(1 - B2q2)5/2 (1 _ B2q2j^ 


lim 

0-^0 



B9(MBq + MBC)\/Bq(BC - Bo) 
1/1 - B202(Ba - B0)^ 


(MBo + NBC)yfeo(BC - Bo) 

/l - b202(i . B2q2j(Bo - B9 J 


2(MBq + NBCJi/bo(BC - Boj/l - b202 


(bo - Be;^ 


nBC 

d(Bo) + / 

'^B(0+t)) 


B0(MBo + NBCJ\/Bq(BC - Ba) 


/l - b20 2(bo - B0 / 


(MBo + NBC)/Bo(BC - BoJ 
/l - b202(i _ B2q2)(Bo - B0 } 


2(MBq + NBC;/Bq(BC - Bcj/l - B^^ 
(Bo - B0p 


d(Bo; - 


2B0(MB0 + NBCJ^(B0 - BC) + B0BC(3M - 2H; + MB^C0|/l - B^^ 


Bti/i - 6^02 


Bth/B0(B0 - BC) 


( 22 ) 
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By performing the integrations in equation (22), the following 
expression is obtained: 


d(v/x) _ rt 

^ /l + k2(l - k2)^^ 


&K'(k)(M + Mk^ + 2Nk) + 


E'(k)(2Mk^ - 2M - Nk - 2Mk^ - Nk^^ (25) 

Consider the rolling case; that is, M = Mp, N = Np, and from 
equations (2) and (4) v/x and are, for 0=0, equal to 0 

and p, respectively. Solving equations (20) and (25) for and Np, 


with N = BCNp = 


2k 


1 + k^= 


Np, gives 


Mp = 


|/l + k2 k2(l + k2)K'(k) + (l - 4k2 + k^)E'(k)^ 


rt[^k^K'(k)2 + k2(l + k2)K' (k)E' (k) + (2 - k2)(l - 2k2)E'(k)^ 


(24) 


Np = 


p(l + k2)^^^[ik2K' (k) - (1 + k2)E'(k)] 

2n[^k^K'(k)2 + k2(l + k2)K* (k)E' (k) + (2 - k2)(l - 2k2)E* (k)^ 


(25) 

For the yawing case M = N = — ^ Nr, and from equations (5) 

1 + k2 

d(v/x) 

and (6) - = -r and = 0. Solving equations (20) and (25) 

^ ae 

simultaneously after making these substitutions yields 


Brk\/l + k2[ik2K’(k) - (l + k2)E'(k)|] 


(26) 


rt[^k^K'(k)2 + k2(l + k2)K'(k)E'(k) + (2 - k2)(l - 2k2)E'(k)^ 


Mr = 


N 
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( 2 ?) 


It is convenient for plotting purposes and in expressing the aero- 
dynamic coefficients to make the following definitions: 


so that , Nj. ' , Mp ' , and Np ' are functions of BC only . The 

variations of these fo\ir parameters with BC are shown in figure 4. 

The velocity potentials for the rolling and yawing motions, com- 
pletely defined hy the velocity potential (eqs. (l6) and (17)) an<i the 
constants given in equations (2l), (24), (25), (26), (27)> and (28), 
may now he written as 




(28) 




0p = px2(h^’0 + Np'C)^0(C - 9) 


(29) 


and 


0r = Brx2(Mr'e + Ni.'C)i/e(C - 9) 


(50) 
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The pressure coefficients for the rolling and yawing motions found 
from equations (29)^ (30)^ (lO)^ and (ll) are 


and 



2p Mp'Cz2 + 5 c2Np’xz - 2CNp'z2 


V 



z) 



2Br Mr'Cz2 + 5c2Nr.'xz - 2CNr'z2 
^ /z(Cx - z) 


(31) 


(52) 


Determination of the Surface-Press\ire-Distribution Expression 


for Constant Lateral Acceleration 

The lateral-acceleration p motion is time dependent and is not 
governed by equation (l) but by the linearized partial-differential 
equation for \insteady supersonic flow; 


2 ^20 d^0 ^ 2V _j_ 1 

8x^ 8y^ 8z^ 8t^ 


( 55 ) 


The boundary condition for the p motion on the tail surface, which is 
approximately in the y = 0 plane, is 


V 



(54) 


The potential function satisfying equations (35) and (54) may be obtained 
from the a potential function given in reference l8 for a wing in the 
xy-plane which was in turn obtained from an analysis by Gardner (ref. 19). 
This potential function is 





( 35 ) 
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where f is the steady-state potential corresponding to a \init yawing 
velocity about the z-axis and X is the steady-state potential corre- 
sponding to a unit angle of sideslip. Thus it is seen that the potential 
function for the p motion may be written in terms of the time-independent 
solutions already obtained. It should be noted that equation (55) dif- 
fers from the formula given in reference l 8 by two negative signs, one 
before the whole expression and one before the second term within the 
brackets. These two sign differences are necessary to account for the 
fact that the boundary conditions for the p and p motions 

V = pV 

V = pvt 

are opposite in sign to those for the a and a motions 

w = -aV 
w = -aVt 


These sign differences on the boundary conditions for the analogous 
motions require that their potentials and pressures be of opposite 
sign. 

The pressure distribution for a time-dependent motion at the 
time t = 0 from the linearized Bernoulli's equation is 



Equation ( 36 ) in conjunction with equation ( 55 ) yields 



(56) 


( 37 ) 
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The pressure coefficient for the sideslip motion may be 
reference 11 as 


2pz 

n , 

BE'(k)\/z(Cx - z) 



obtained from 


(38) 


The reciprocal of the complete elliptic integral E'(k) appearing in 
equation ( 38 ) is plotted in figure 4 for various values of leading-edge 
sweep parameter BC . 


The expression for X, the only quantity in equation (3?) as yet 
undetermined, can be obtained by substituting equation ( 38 ) with 3=1 
into 




r 

AP(x2^,z) 

l.E. 

__ 1 _ 


dxi 


-13=1 


(39) 


and carrying out the indicated integration. This process gives 


X 


VV2(l - ^1 - B^^)</z(Cx - z) 
BCE'(k) 


(40) 


The pressure coefficient for the 3 motion is now completely defined by 
equations (52), ( 58 )^ (57); (^O). 


Force and Moment Coefficients for 
Zero-End-Plate Triangular Tails 

With the pressure distributions known (eqs. (51), (32), (37); and 
( 58 )) the total side force^ yawing moment about the apex, and rolling 
moment about the root chord can be determined for the various motions 
by the formulas 


Side force = 



(41) 
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Yawing moment 



— X dz dx 

q 


(h2) 


Rolling moment 



( 45 ) 


The force and moment coefficients and the stability derivatives may be 
readily obtained from these quantities for zero-end-plate triangular 
tails and are given in table I. 


Force and Moment Coefficients for Complete -End-Plate 

Triangular Tails 

For the sideslip, yawing, and constant- lateral-acceleration motions 
the pressures acting on the vertical tail in the presence of a complete 
end plate are the same as the pressures acting on one-half of a symmet- 
rical wing for the angle-of-attack, pitching, and constant-vertical- 
acceleration motions, respectively. Each half of this symmetrical wing 
should have the same plan form as the vertical tall under consideration. 
For the sideslip motion the pressiires on the vertical and horizontal 
tails with a partial end plate have been reported in reference 9* 

Attention is now directed to the effect of end plates on the sta- 
bility derivatives. For the case of the complete end plate the following 
vertical- tail derivatives can be obtained from symmetrical- wing results 
(ref. 20): and Cn^. The transformations 

needed to change the symmetrical-wing derivatives into these vertical- 
tall derivatives are: 


Expression for Cyp = -(^Expression for with A replaced by 2Ay) 

Expression for Cn„ = X (Expression for Cm with A replaced by 2AY") 
p 5Av \ ^ / 


Expression for Cy = X (^Expression for Cl with A replaced by 2Ayj 
^ 5At' V q ^ 


Expression for Cn = 5 X ^Expression for Cm with A replaced by 2Ay^ 

^ 9Av 
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Expression for Cy^ = X ^Expression for with A replaced by 2Ay^ 


Expression for Cn^ = — X ^Expression for Cm^ with A replaced by 2Ay^ 


It should be noted that the above coefficients obtainable from wing 
results do not include the rolling derivatives ^Ij.} ^^3’ 

Before these aerodynamic coefficients can be evaluated, the rolling 
moment induced on the end plate must be determined. This moment for 
the 3 motion has been derived in reference 9 the induced moment 
for the other two motions are not yet known. 


RESULTS AHD DISCUSSION 
General 


Table I contains the formulas for all the p, r, 3 , and 3 stability 
derivatives for the tail shown in figure 1. Table II presents these 
derivatives for the half-delta tail (the special case of zero trailing- 
edge sweep, AyB = 2BC). The variations of these half-delta derivatives 
with the parameter BC have been plotted in figures 5 "to 10. It is 
evident from tables I and II that the expressions for the ^ derivatives 
have been separated into two component parts, each part being multiplied 
by a different function of the Mach number parameter B. Each of the 
components, excluding this factor, is a function of BC. Figures 6, 7 ^ 
and 8 show the variation of the two parts with BC for each of the 
3 stability derivatives. Once the Mach number, and hence B, has been 
specified, the two parts may be combined and the total derivative deter- 
mined for any given aspect ratio or leading-edge sweep . For the reader ' s 
convenience the variations of the 3 stability derivatives with Mach num- 
ber for a n\imber of aspect ratios of the half -delta tail have been plotted 
in figures 11, 12, and 15 . Mach nixmber variations of the other half- 
delta stability derivatives may be obtained by inspection from their 
variation with BC . 

Side-force and yawing-moment coefficients (as obtained from ref. 20 ) 
of the 3, and r motions for the complete -end-plate vertical tails have 
been presented in table III. As in the case for the zero end plate, the 
complete-end-plate stability derivatives for the half-delta tail are 
presented in a separate table (table IV). The quantities 1 /e'(BC) and 
G(BC) appearing in the expressions for the complete-end-plate stability 
derivatives have been plotted in figure li)-. 


NACA TN 3240 


23 


It should he realized that the force and moments given by equa- 
tions (4l), (4-2), and (4-3) and the stability derivatives in tables Ij 
II, III, and IV are for axis locations at the tail apex and root chord. 
Formulas for the transfer of force and moment coefficients to a body 
system of axes with the origin displaced a distance Xq (positive 
forward) from the tail apex and a distance Zq (positive downward) 
from the tail root chord (see fig. 3(h)) are presented in table V. 


End-Plate Effect for Sideslip, Yawing, and 
Constant-Lateral- Acceleration Motions 


Figures I 5 to I 8 have been prepared to show, for the half -delta 
tail, the effect of a complete end plate on the side-force and yawing- 
moment aerodynamic coefficients. It is evident from the large magnitude 
of the differences between the two limiting cases of a complete end 
plate and no end plate shown in these figures that a reasonable estima- 
tion of the partial-end-plate effects would be highly desirable. End- 
plate effects for various sizes of horizontal tails for the sideslip 
motion have been evaluated exactly in reference 9 hut are not presented 
herein. The following approximate formulas for the side- force aind yawing- 
moment stability derivatives based on these partial-end-plate results 
are suggested: , 


(<^Yr)p ■ (%)c 




where 


('^nr)p - ('^'v)c ■ 1 

0 

1 

(%)p ' (%)c - [ 

(%)c ' (%)o 

1 

0 

a 

II 



(%)c-(%)p 
' ’ - (%)o 


(Wt) 
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It has been mentioned previously that the effect of end plates on 
the rolling-moment coefficients has only been evaluated to date for the 
sideslip motion (ref. 9)* By using the results of reference 9, figure 19 
has been prepared to shov the variation with BG of the stability deri- 
vative Cip for both a half-delta vertical tall mounted on a con^jlete 

end plate and an isolated half-delta vertical tail. Also plotted in 
this figure are the vertical-tail and end-plate contributions to the 
Cjp derivative of the complete-end-plate — vertical- tail combination. 

Partial-end-plate effects on the rolling-moment coefficients of a 
vertical tail perfonnlng ^ and r motions cannot be approximated by equa- 
tions similar to equation (^A-) since the contribution of the complete 
end plate to the total rolling moment for these motions is unknown. 


End-Plate Effects on Rolling Motion 

In the analysis of end-plate effects, conqjlete-end-plate stability 
derivatives were evaluated by using wing results. This treatment was 
possible because in the sideslip, yawing, and constant- lateral- 
acceleration motions the complete end plate acts only to uphold the 
loading in the same manner as one half of a wing does on the adjacent 
half. The complete end plate on a rolling vertical tail with, of course, 
the end plate rolling with the vertical tail causes a sidewash in the 
plane of the vertical tail which can result in large induced loads, 
depending on the size of the vertical tail relative to the end plate. 

This is an end-plate effect very different in nat\are from the one experi- 
enced in the sideslip, yawing, and constant- lateral-acceleration motions 
where the end plate does not cause a sidewash in the plane of the vertical 
tail. Clearly then the complete-end -plate boundary can only be esta- 
blished by solving the difficult nonplanar problem. However, some idea 
of the end-plate effects may be obtained by considering the results 
shown in reference 5 slender nonplanar tails . Reference 5 shows 
that for ratios of vertical-tail span to end-plate span greater than 
0.75 the vertical tail does not experience any large changes in loading 
due to the end plate. This fact suggests the possibility that the Mach 
number and aspect- ratio effects on the isolated rolling tail (fig. 9) 
might be applied to the slender, nonplanar, vertical- tail loadings of 
configurations with ratios of the vertical- tail span to end-plate span 
greater than O.75 to yield good estimates (from a theoretical viewpoint) 
of the side force and yawing moment. 

Aspect-ratio and Mach number effects on the rolling isolated tall 
are apparent from figure 9* The same effects for a rolling wing, which 
may be thought of as an end plate with a zero- span vertical tall, are 
illustrated in reference 13. In both cases the magnitude of the forces 
can be predicted within 10 percent by slender theory up to leading-edge 
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sweep peLTajueters BC of 0.5* is not unreasonable to expect that 
the same agreement would exist for nonplanar tail arrangements as long 
as the sweep of the leading edges of the panels of the configuration 
satisfy the condition 

(B X Tangent of apex angle of panel) <0-5 


Lateral Oscillatory Motion 

It has long been the practice to estimate the aerodyna m ic daji 5 )ing 
of longitudinal low- reduced-frequency oscillations of lifting wings from 
the damping in pitch and damping of the vertical-acceleration motion. 

An analogous approximation may obviously be used for isolated vertical 
tall surfaces oscillating laterally in yaw; that is, the lateral damping 
can be estimated from the results obtained for the yawing and constant- 
lateral-acceleration motions. This approximation which represents the 
first-order frequency terms is given by Cn^. - Cn^- 


Figure 20 has been prepared to show the variation with Mach number 
of this damping (^Cnj. - Cn^) for half-delta tails of aspect ratios 1.0, 

1.5, and 2.0 with a conq)lete end plate and with no end plate. The yawing 
axis for these examples is located at the tail apex. In order to illus- 
trate the relative magnitudes of the two tenns, the Cnp contribution 

is also plotted in figure 20. It should be kept in mind in comparing 
the no-end-plate damping with the complete-end-plate damping that Cn^, 


is always greater (more negative) for a vertical tail with a complete end 
plate than with no end plate and that for any given aspect ratio Cn^. 

will decrease (become less negative) as the Mach number is increased. 
(|^Note that negative values of Cn^. ~ indicate positive damping.^ 


For the Ay = 1.0 tail shown in figure 20(a), the total damping 

of the complete-end-plate vertical tail is greater than that of the no- 
end-plate vertical tail. This occurs even though the contri- 

bution for the complete -end -plate tail is positive and therefore 
detracts from the total damping, Cn^ + whereas the con- 

tribution of the zero-end-plate Ay = 1.0 tail is slightly negative 
over the Mach number range for which the theory is valid and hence adds 
to the damping. As the aspect ratio is increased from 1.0 to 1.5- 
(fig. 20 (b)), the damping of the complete-end-plate vertical tail 
decreases at a more rapid rate than the danq)ing of the no-end-plate tail 
to the extent that the total damping for the corap lete-end-plate Ay = I.5 
vertical tail is now slightly less than the no-end-plate tail. This 
change may be attributed to the fact that the contribution to the 
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complete-end-plate tail dancing increases more rapidly (becoming more 
positive) than the ("^np^ con^onent of the no-end-plate tail. Figure 20(c) 


for the Ay = 2.0 vertical tail shows that the damping of the complete- 
end-plate vertical tail continues to decrease more rapidly than for the 
no-end-plate vertical tail so that the complete-end-plate tail damping 
becomes considerably less than that of the no-end-plate tail. 


The damping derivatives plotted in figure 20, as pointed out, are 
for an axis location at the tail apex. In order to depict the effect 
of moving the yawing axis forward from this point, dandling derivatives 
have been computed for the aspect- ratio- 1. 5 tail with a yawing axis 
located 1 chord ahead of the tall apex (fig. 2l). These contputations 
have been made with the aid of the transfer formulas in table V. A 
comparison of figures 20(b) and 21 shows that the damping qualities of 
the complete-end-plate vertical tail relative to the zero-end-plate 
vertical tail were considerably improved by moving the yawing axis 
forward. This improvement can be acco\mted for by the appearance of 
the p derivatives in the transfer formula for Crij,. 


C0NCIUDING REMARKS 


Pressure-distribution expressions and stability derivatives have 
been derived for zero- end-plate triangular vertical tails performing 
yawing, rolling, and constant- lateral- acceleration motions by a method 
for solving supersonic- conical- flow boundary- value problems. In addi- 
tion, by using the yawing and constant- lateral-acceleration results, 
the damping of a vertical tail oscillating laterally in yaw has been 
approximated to the first order in frequency. End-plate effects have 
been discussed 6ind suggestions made to aid in their evaluation. In 
this connection the complete-end-plate and no-end-plate stability deri- 
vatives for the sideslip motion obtained from other sources have been 
considered. 

The pressure-distribution expressions and stability derivatives 
contained in this report are valid for a range of Mach numbers for which 
the leading edge is subsonic and the trailing edge supersonic. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 10, 195^. 
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APPENDIX A 

DEVELOPMENT OF EqjJATIONS RELATING THE v- VELOCITY TO THE DISTRIBUTION 
FUNCTION f(a) IN THE y = 0 PLANE 


Equation (8) gives the expression for the velocity potential every- 
where in space resulting from a distribution of doublets in the xz-plane 
with the strength of each doublet in this distribution being governed by 
the doublet- strength function A(x,z). The derivative of this velocity 
potential with respect to any one of the coordinates x, y, or z will 
give the pertiirbation velocity in that direction. Of primary interest 
herein is the v-velocity, or the y-derivative of this potential, that 
is, 


v(x,y,z) 




(Al) 


for points on the xz-plane. Brown and Adams in reference 15 have con- 
structed the velocity potential in space of a distribution of doublets 
by use of the following procedure. First, by using equations (8) and 
(l8), the potential of a line of doublets in the xz-plane at an angle 

tan‘*^(J to the x-axis is detemilned. This potential is given by 


0 = 


B2y(x - _ 

(1 - 


j \ 2B% /x^ - B^(y2 + z2) 

- V ^ (1 _ B2a2)^ 


(A2) 


where 


X - B^az 

01 - B2ff2Yx2 _ B2(y2 + z2) 


The velocity potential of a distribution of line doublets in the xz-plane, 
on the vertical tail, with strengths governed by the distribution func- 
tion f(a) may then be written as 


0 



f(a)0L da 


where 


tan"^ C = e, the apex angle 


of the vertical tail. 


(A5) 
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Substituting equation (A2) into equation (A5) and differentiating 
with respect to y yields the following equation for the v-velocity 
as py/x approaches zero: 


y 

X 


lim ■<’ 

Ti-^0 




'0 


Bf(g)V^l - B^^(l - B^oe 
(1 - B^a^)^(Ba - B0)^ 


3Bf(g)(l - B^oe)coth-l| 
(1 - B2g2)5/2 


2Bf(g)\/l - B2e2 

(1 - B2g2)2 


d(Bg) + 


rBC 

^B(9+t]) 


Bf(g)|/l - b202(i _ B^ge)^ 
(1 - B^g 2 )^(Bg - B9 


5Bf(g)(l - B2ge)coth-l| 2Bf(g)|/i - b202 


(1 - B2g2) 


5/2 


(1 - B2g2)‘^ _ 


d(Bg) - 


2f( e)/l - b2&2 


(A4) 


The singularity vhich occurs in the term of equation (A2) when 

72 - 1 

y is set equal to zero has been accounted for in equation (Ak) (see 
ref. 13 )• 

By taking the first and second derivatives of equation (a 4) with 
respect to 0^ two other useful relations sire obtained. They are given 
in the appendix of reference 15 as 


de 



3B^af (a)coth*‘^| 

(1 - B2a2)5/2 


B^(3Bg + 2B8 + BeB^q^)f(g) ^ B8B2f(q) 

/l - B^^(l - /l - B^2(Bo - B9)^ 


B2f(o) 


2B2f(a)/l - b202 


1/1 - b 20 2 (i _ B 2tj2)(Ba - Be) (Bo - B6 )5 


rBC 

d(B0) + / 

-D / 


B(e+ii) 


3B^qf ( a ) eoth~ 
(1 - b2c2)5/2 


b2(3Bo -t- 2B6 + BeB2g2)f(o) ^ B9B2f(g) B^f(a) 

/l - b202(i _ B2g2)^ /l - b262(bo - B0 - B%^(l - B2g2)(Bc - Be) 


2S2f(g)/i - b2b2 

(Be - B6)5 


d(Bo) - 


2B2ef(e) u/i - B^^f'(e) 

Tl/l - b 262 n 


(A5) 
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and 


d^(v/x) 



(Ba - B0)^ 


d(Bff) 


+ 


6/l - b202 


■'BC 


B5f(a) 

B(e+T)) (Ba - B0) 


^ d(Ba) - 


l/l - b202 


6f”(e) 

n 



(a6) 


The factor multiplying the 

as It appears in reference 
in equation (a6). 


f"(0) 


d^(v/x) 

term of the expression for 

ae2 


13 is slightly in error and has been corrected 


Considering equations ( 5 ) and ( 7 ), it is evident that equation (a 6) 
must be zero for both the rolling and yawing cases. This equation has 
already been satisfied by f(0) (eq. ( 17 ))^ since equation (a 6) is in 
essence the integral equation which was inverted to obtain the general 
pressure expression from which f(0) was derived (see ref. 15 )* 
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APPENDIX B 


INTEGRATIONS TO OBTAIN v/x 


Z = i 

X B 


The expression for v/x is, for 0 approaching 0 (see eq,. (l9))> 

Q 

-5(MBct + NBC)/Ba(BC - Bcr)tanh"^/l - B^cr^ 


(1 - B2a2)5/2 


© 


2(MBff + NBC)/Ba(BC - Ba) 



(l - B2g2)^ 


f 

oB(0. 

lim < 

lim 

i / 

0 -^0 

q-^O 

V. 


pBC 

(MBct 4- NBC 

^B(0+ti) 


d(Ba) + 


(MBq + NBC)\/Bg(BC - BaXl - B^o0)^^^ B%%(Bg) 

(l - B2g2)2(Ba - B0 


(1 - B2g2}^(Bg - B0)^ 


© 


2(MB0 + NBC)\/b 0(BC - B0)|/l - b2©2 




(BIJ 


This expression has been broken into parts as indicated by the circJ^d 
numbers with the third part being broken into two additional parts 
and because of the singularity in the integrand. Since and (3^ 

are elementary integrations similar to those found in most Integral 
tables (see ref. l6), only (^ is dealt with in detail. 
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Performing the integrations and and combining the results yields 

j 2(MB9 + NJ>/B9 (bC - B9 J 
B I Bti 

Jt |M(-7EeBC - 2 + 10B9 - BC) + N(4 - 7BC + 4B9 - B9BC J ^ 

8/l - BC(l - B9) 


rt[M(-2 - 1OB0 + BC - 7B0BC) + N(-l4- + 4B9 - 7BC + B0BCjJ 

8/L + BC(l + Be) 


(B2) 


where N = NBC. The first term of expression (B2) exactly cancels 
and the total of ® 


Jt 

8B 


M(2 + BC) - N(4 - 7BC) ^ M(-2 + BC) - N(4 + 7BC) 


|/l - BC 


/l + BC 


(B5) 


The following two integrals comprising remain to be evaluated: 


-3 MBa/Ba(BC - Ba)tanh“®/i - B^a^dCBc) 


B 


(1 - B^a^) 


2 0.^2 


(Bk) 


_3 N^a(BC - Ba)tanh-l/l - B^a^(Ba) 

) (1 - b2c2)5/2 


B 


(B5) 


It should he mentioned at this point that the integrands of expres- 
sions (b 4) and (B5) are finite and continuous over the interval 0 to 
BC and therefore must yield a finite quantity when integrated. 
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The integration of expression (b 4) hy parts gives 


-m/Bct(BC - Ba)tanh-ll/l - S2a2 


B(1 - B2a2)^/^ 


BC 


B 


0 


(BC - 2BaJtanh"ll/l - 

2(l - B2a2)^^^/Ba(BC - Ba) 


|/Ba(BC - Bo)d(Ba) 


Bo(l - b2o2) _ 


'0 


(b6) 


Integration of expression (B5) hy parts gives 

iBC 


-N^a(BC - Bajtanh-^/l - B 2 ff 2 


_ BB0(1 - b2o2) 


3/2 


3C - 


->0 


KBC tanh“^/l - B^a^d(Ba) 

0 2BBa(l - B2o2)^/^\/fea(BC - Bo' 


n/bo(BC - Ba)d(Bo) 
0 BB2 o 2 (i - B 2a2)^ 


(BT) 


Combining expressions (b 6) and (B7) results in 




-(MBo + N;\/bo(BC - Bo)tanh"li/l - 


BBo(l - b2o2)^/^ 


BC 


(MBo + N)^o(BC - Bo}d(Bo) 
3 BB2a2(i _ s2a2;^ 


[^o(BC - 2Bo) - NBC^tanh'^/l - B2a2d(Bo) 




2BBo(l - B2a2j^^^y^a(BC - Bo) 


(b8) 


The first term of equation (b 8) ^ when evaluated at the limits^ is either 
zero or infinity. The integrand of as was noted is finite over the 
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whole interval; therefore^ infinities introduced as a result of parts 
integrations must, in the end, cancel themselves. 

The second term of equation (b 8) is an elementary integration which 
when evaluated (with infinities neglected) yields 


pM + 4N - 3MBC - 3KBC 
8 b\ /l - BC 


-4N - 3MBC + 2M + 3MBC 
/l + BC / 


(B9) 


It is now convenient in integrating the third term in equation (b 8) 
to introduce the variable substitution 


Ba 


+ k 
1 + 


so that BC 


and k are related by 


(BIO) 


BC 


2k 


1 + k^ 


(Bll) 


The third term in equation (b 8) when transformed by equation (BIO) may 
be written in the form 


^ li (B12) 

Byi + k2(l - k2) 

where 

^-k 1 - r 

l2 = - / k(M + Nk5)F(U<i^ 

^-k 

1 - 
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I, = Yam 


-k 


h 


r^ 3Nk?^ 

L 


F(Ud^ 


r^ Nk2(l _ k2j . ^ ^ 

I6 = / ^ F(UdC 

^-k (1 - ^2)(k2 . ^2) 


Irj - 


Nk(l - k4^2^^ 

J-k (1 - ^2)(k2 _ ^2) 


F(UdC 


ta^-lfAUjO! 

\ 1 + k^ 

7(1 - ?2)(k2 - ?2) 


The integrals and ly are elementary and may be determined by 

integration by parts . If the multiplicative factor before the summation 
sign in equation (B12) is neglected until all the components are totaled, 
these three components become 


ii + 15 + ij 


-Mk(l + k4jjt mc2(i + 

1 - k2 1 - k^ 


Consider the integration required for I 2 , that is. 


(BI 3 ; 


tanh" 




-k 


• i/(i - 

1 -t- k^ 

]/{l - ^2j(k2 . ^2) 


d^ 


(Bl4) 
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Let ^ = k sin 9; then expression (B14^ becomes 


nrt/2 

tanh"^ 

\J {1 - k2)(l - k^sin^O ) 

_ 1 + k2sin 0 _ 

J -n/2 

1/1 - k2sin20 


d0 


(B15J 


It can be shown that 


tanh" 


\j(l - k^)(l - k^sin^e) 


1 + k^sin 0 _ 


= tanh 


-1 


]/l - k^ 

/l - k^sin^e 


- tanh' 


_1 \ll - k^sin 9 
/l - k^sin^e 

(Bl6) 


This fact allows expression (BI5) to be rewritten as 


f 


, rt /2 


tanh" 


- rt /2 


/l - k2 

\Jl - k^sin^y 

^1 - k^sin^O 


d0 


nJt 




/2 


tanh' 


/2 


_n/|/l - k^sin 0' 
- k2sin2 0y 
/l - k^sin^O 


d0 


(BI7) 


The integrand of the second integral is an odd function; therefore, 
the integration between -jt/2 and rt/2 of this function is zero. Since 
the integrand of the first integral is an even function in 0, this 
integration may be expressed as 



tanh“^| 


\/l - k^ \ 
Yl - k^siii^/ 


d0 


(BI8) 


'0 


\Jl - k2sin20 
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After the inverse hyperbolic tangent is replaced by its logarithmic 
equivalent and the additional variable transformation 


sin^v = 


1 - k*= 


1 - k^sin^Q 

is introduced, expression (b 18 ) becomes 


(B19) 


prt/2 

sin”^/l-k2 


loge 


1 + sin V dv 

^ /sln^v - (1 - k2) 


It is convenient to let 


A = sin"^ ^ - k^ 


(B20) 


then 


sin^A = 1 - k2 


Making these substitutions in expression (B20) gives 



1 + sin V 
1 - sin V 


dv 



sin^A 


(B21) 


which is exactly in the form of the fourth integration formula of 
table 355 ill reference I 7 . This formula gives the value of expres- 
sion (B21) as rtK'(k). The integration of I 2 may now be expressed as 

I 2 = -k(M + m^)rtK'(k) (B22) 

Using the same integrating procedure for and Ig as just 

outlined for I 2 and the integration formulas in tables 555 556 

of reference I 7 leads to 

+ k2) + 3 Nk^jt [1 + K'(k) - E' 

1 - k2 




•NACA TN 3240 


57 


% + I6 = jtN(l - k2)k2n |-(l-k2j, - N(l - k^jk^nK' (kj - 


Nk2rt |l + k' (k) - E' (k£] 
1 - k2 


where II |-(l-k2), ^-k2,i | is a complete elliptic integral of the 
third kind with modulus |/l - k^ and parameter -(l - k^j. 


Summing all the various parts contributing to the third term in 
equation (b 8), including the common factor, gives the following 
expression: 


-1 


B(l - k2)|/l + k2 


rtMk^ + nHk^ + 


jtK'(k) [ik^M + Nk2(l + k2)]] 


1 - k2 


«E’(k)f2k2N + Mk(l + k2)l - ^ , o 
!= =!■ + Nk2n(l - k2)n 


(B25J 


1 - k^= 


The addition of expression (B23) to expression (B9) completely 
evaluates (^ • Expression (B5) gives the evaluation of CD^ 

(^. Before writing the total integration, the sum of expressions (B25 )j 
(B 9), and (B5)^ it is desirable to combine expressions (B3) and (B9)^ 
which are fxinctions of BC, and transform them by equation (Bll) to 
functions of k. This procedure yields 


jt(Nk2 + Mk5) 


(B24) 


b/i + k2(l - k2) 
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The total integration may now he written in terms of the param- 
eter k as 


b/i + k2(l - k2)‘ 


K'(k) [ik^M + Nk2(l + ]^^)] - E’(k) [ik^N + Mk(l + k2)J + 


Nk2(l - k2)^n 


(B25) 


By use of the process commonly known as interchanging the amplitude and 
parameter (see pp. 155 to l4l of ref. 21) the elliptic integral of the 
third kind appearing in equation (B 25 ) is found to be equivalent to 

This operation permits the expression for v/x (eq. (B25)) to 
k^ 

assume the form given in equation (20) in the body of the report. 
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TABLE I.- LATERAL STABILITY DERIVATIVES FOR ZERO-END-PLATE 
TRIANGULAR VERTICAL TAILS 


Derivative 

Formula 

»(%)o 

^ VAyB\ 4 2AvB / 


^AyBj^^Mp’ , 5Np' ^ 5BC»p'\ 
16 Kb\ 2 ^2 * AyB J 

(S)o 


^ V52BC l6AvBy ^ \^32BC 4AyB Qp,^2-q2 J 


(%)o 

^ 1/a^ 4 4 2AyB / 

o 

«AvB ^/5Mr ' 5Nr ' 5BCNr 

16 1(AvB\^ 2 ^ 2 ^ A^B y 

(^nr)o/® 


^ (32BC leAys) (32BC ■" l+AyB ■" SAy^^) 


K%')o 

-It \/l - \ll - B^C^ 

E'(k)\/BC 

»(C.p)o 

-it^/A^V^ - \/l - B^C^ 
2E'(k)v^ 

(Cnp)o 

- )ll - b2c2/5AvB ^ 
5B2c2E'(k)|/A^ ^2 / 


k2 
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TABLE I.- LATERAL STABILITY DERIVATIVES F® ZERO- END-PLATE 
TRIANGULAR VERTICAL TAHS - Concluded 
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TABLE II.- LATERAL STABILITY DERIVATIVES FOR ZERO- END-PLATE 
HALF-DELTA VERTICAL TAILS 


Derivative 


Formula 


“(%)o 


+ 2 Hp') 


(S)c 


-« |(Mp' + 2 Np'; 






(CYr)o 


5|2(Mr' + 2Hr') 


(‘'“r)o 


^(Mr' + aHr') 

o 




§ M/ . Mr') 


B(CYp)o 


- l/l - 

E'(k) ^ 


(S)o 


— \L(l - l/T^ B^C^) 

3BCE'(kr 


B(Clp)o 




2E'(k) 




. 2 Nr' -4 

3B2c2E'(k; 


BCjt 

2B^ 


Mr' + 2Hr' - 2 


- \/l~- 
B^C^E'Ck) 




('^2p)o 


kB 

4 


B^C^E'Ck) 


I Mr' + 3 Nr' - 3 


\4(i - 1^- B%:^) 


B^C^E'Ck) 


-BCrt 

4 


BCn 

4B^ 


2Nr' 


^ ^ 1^1,- 

4- otn2r*2T? ' ^ i, » 


2Hr' + - Mr' - 2 


2B2c2E'(kj _ 

|4(l - \[T~- B^^) 

B^C^E'Ckj 
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TABLE III.- LATERAL STABILITr DERIVATIVES FOR COMPLETE-EHD-PLATE 
TRIANGULAR VERTICAL TAH^ 

[Formulas obtained from reference 2(^ 


Derivative 


Formula 






» .1.-% . biST^) 

E'(BC)(1 + 


(^p)c 


4(1 - R)^ (7 


5E’(BC)(1 - r2) 


,5/2 


1^2 + R^)(| + Bln-^) + R(4 - r2)/i - r2 


fYr)c 


4g(BC)( 1 - ^ sln-^) + R(5 - 2 r2)|/i - r2 

3(1 - r2)5/2 ^ ^ 


(C^r)c 


-BG(BC ^ 1 - R p ^ 5)(| + + R(4 r‘^ - 12r 2 + 23)|/l - R^ 

3BC(1 - r2)'^/2 L 


("^p)c 


4(1 - R)^ 
5(1 - r2)^/^ 




b2e'(BC) L-V 

G(BC)[5(^| + sin-lR) + R(5 - 2 r2)/i - r2 + 

1^2 + r 2)(| + sln-lR) + R(4 - r2)i/i - r2 




-(1-R)5 I 

^(1 - e’(bc)g(bc)] r 

3(1 - r2)'^'^Sc' 

, BE’(BC) [ 


12r2 + 23) /l - r2 + B ^ 3(^1 + sln-^) 


^-G(BC)(2r2.3)|. 




R|/l - r2 


^ ~ ^ - G(BCK^R*^ - 1®2 + 23 ) 


E'(BC) 
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TABLE IV.- LATERAL STABILITY DERIVATIVES FOR COMPLETE- END-PLATE 

HALF-DELTA VERTICAL TAIIS* 



1 - 

(1 - 2b2c2)e'(BC) + b2c2k'(BC) 


Derivative 

Formula 

K%)o 

-2kBC 

e'(bc) 

(%)c 

kn 

3e'(bc; 

(^Yr)c 

2jtG(BCj 

(^nr)c 

-5«bg(bc; 

2BC 

(%)c 

^ PcfEC ) + ^ 

p 

1 rfvr ) 4- 

b2Q^ e’(bc;_ 

T Cjl 1 

_ 3E (bc; J 

(S)c 

TcfEC) 1 ^ 

1 ^®rGCBCj+ ^-1 

2BBC [_ E (BC )_ 

2BC j_ 5E ' (BC J 


*(cip)c 


can “be obtained from reference 9* 
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TABLE V.- TRANSFER FORMJLAS 


Stability derivatives in a 
body system of axes with 
origin at tail apex 

Formulas for transfer to a body system of 
axes with origin displaced distances Xq 
( positive forward) and Zq (positive 
downward) from tail apex 

CYp 



"“P - ^ % 


^ % 

% 

% 

Cnp 

""p - ^ % 


c, + -^ Cy 
‘P by P 

Cy, 

^ % 

Cn 

xir 

V - ^(S " "^0 

"2r 

Cl + -^ Cy 2 do - -2 Cy 

r t)^ r by P by by 3 

% 


Cnp 

^ % 

^2p 

"'p " ^ % 



Figure 1 .- Sketch of the vertical tail showing the axes system used in L-8562I4. 
the analysis and the positive directions of the forces, moments, 
velocities, and angles. Dashed lines on vertical tail indicate that 
the trailing edge may he swept forward or backward. 
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End plate 


Vertical tail 


Figure 2.- Half-delta tail mounted on a complete end plate. L -85625 
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(a) Positive direction of forces and moments in the stability axes 

system (a = 0^). 


z 



(b) Axes system to vhich force and moment derivatives may be transfenred 

by use of table V. 

Figure 5«- Systems of axes. 
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Figure 4.- Variation of the parameters Np'BC, Mj.', Np', Mp', and 

l/E'(k) with BC. These curves are of use in computing the no-end- 
plate stability derivatives given in tables I and II. 




Figure 5.- Variation of the zero-end-plate half-delta vertical- tail 
stability derivatives (pYr^Q C^^r)o stability- 
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derivative parameter 
parameter BC . 
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with the leading-edge sweep 
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Figure 7 .- Variation with BC of the quantities 


=«[(S)o]; 


and 


® l(S)o ]2 for the half-delta tail. Note that the stability 
derivative is given by the sum of ^ and 
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Figure 8.- Variation with BC of the quantities 


for the half-delta tail. Note that the stability deriva- 


tive 




0 


is given by the sum of 
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Figiire 10.- The variation of -B^CYp^Q^ and (^C^p^Q with 

BC for the half-delta tail. 
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Figure 11.- The variation of 

aspect ratios of a half-delta tail 
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Figure 12.- The variation of with Mach number for a ninnber of 

aspect ratios of a half-delta tail. 
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Figure 15.- The variation of 


(^*3)0 

aspect ratios of a half-delta tail 


vith Mach number for a number of 
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Figure l4.- Variation of the parameters 1 /e'(BC) and G(BC) with BC. 
These curves are of use in computing the con^lete-end-plate stability 
derivatives given in tables III and IV. 
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Figure l6.- Complete-end-plate and no-end-plate boundaries of the 

quantities [(%)ll* Note that the derivative Cy* 

is given by the s\m of [^*"^13)12* 
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Figure 17*- Complete-end-plate and no-end-plate boundaries of the 


quantities B 




is given "by the sum of 
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Note that the derivative 
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Figure l8.- Con^ilete-end-plate and no-end-plate boundaries of the 

stability derivative Cn^ and the stability-derivative parameter -BCyp- 
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Figure 19 •- Variation of and the horizontal- and 

vertical- tail contributions to with BC. 
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(a) Ay = 1.0. 

Figure 20.- The variation with Mach number of the first-order approxima- 
tion to the damping of the lateral oscillation in yaw for complete- 
and no-end-plate half-delta tails with aspect ratios of 1.0, 1*5> 
and 2.0. The yawing axis is located at the tail apex. 


OtS^ Ml VOVM 


NACA TN 32^+0 


67 




(c) Ay = 2.0. 
Figure 20.- Concluded. 
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Figui'e 21.- Th6 variation with Mach nuiriber of the first-order approxima- 
tion to the damping of the lateral oscillation in yaw for a complete- 
and a no-end-plate half-delta tail of aspect ratio 1-5 . The yawing 
axis is located 1 chord ahead of the tail apex. 
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